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Abstract 

Stress rupture failure of Carbon Composite Overwrapped Pressure Vessels (COPVs) is of serious 
concern to Science Mission and Constellation programs since there are a number of COPVs on board 
space vehicles with stored gasses under high pressure for long durations of time. It has become customary 
to establish the reliability of these vessels using the so called classic models. The classical models are 
based on Weibull statistics fitted to observed stress rupture data. These stochastic models cannot account 
for any additional damage due to the complex pressure-time histories characteristic of COPVs being 
supplied for NASA missions. In particular, it is suspected that the effects of proof test could significantly 
reduce the stress rupture lifetime of COPVs. The focus of this paper is to present an analytical appraisal 
of a model that incorporates damage due to proof test. The model examined in the current paper is based 
on physical mechanisms such as micromechanics based load sharing concepts coupled with creep rupture 
and Weibull statistics. For example, the classic model cannot accommodate for damage due to proof 
testing which every flight vessel undergoes. The paper compares current model to the classic model with 
a number of examples. In addition, several applications of the model to current ISS and Constellation 
program issues are also examined. 


Nomenclature 

a Weibull shape parameter for composite strength 

P Weibull exponent and shape parameter for composite lifetime under fixed fiber stress 

0 creep exponent for matrix 

p power law breakdown exponent reflecting sensitivity of lifetime to fiber stress 

ct fiber stress 

o p fiber stress at proof pressure level 

CT re f characteristic fiber stress level 

a v composite reference strength (same as cr ref ) 

C, Weibull shape parameter for fiber strength 

F{t) cumulative distribution function for composite failure 
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A 

k ’ same as k 

A 

k critical fiber break cluster size for catastrophic failure 

k p fiber break cluster size after proof 

R(t ) vessel reliability over time 

t time in hours 

t c characteristic time for creep to occur 

t c characteristic time for matrix creep 

t p hold time at proof pressure 

t P} b proof benefit time scale 

t Tof characteristic time scale for fiber failure 

Introduction 

Carbon composite overwrapped pressure vessels (COPVs) are widely used for storing gasses under 
high pressure in a wide range of applications including onboard spacecraft such as the Space Shuttle and 
the International Space Station (ISS). The use of COPVs is currently being planned for the Crew 
Exploration and launch vehicles as well. The principal advantage of the carbon COPV technology over 
conventional all metallic storage devices is the substantial weight savings. The National Aeronautics and 
Space Administration (NASA) has been supporting the development of this technology since the early 
1970s (Ref. 1) with an interest in safe application of these components to reduce mass to orbit. NASA 
Johnson Space Flight Center White Sands Test Facility (WSTF) in collaboration with NASA Jet 
Propulsion Laboratory (JPL) has been testing components in support of this objective since the 1980s and 
WSTF has been involved in test development and analysis to address effects of impact, propellant and 
cryogenic fluids exposure on Kevlar (DuPont)/epoxy and carbon/epoxy composite overwrapped pressure 
vessels (Ref. 2). 

Space flight COPVs are generally made of a thin metallic liner overwrapped with a high strength 
fiber/epoxy composite. This design allows for significant weight savings in contrast to an all metal 
pressure vessel design and these vessels are widely used for spacecraft and commercial applications 
because of this advantage. Mechanisms associated with the failure of the COPV system must be 
considered for safe operation from manufacturing to disposal. The fibrous composite overwrap carries a 
significant portion of the pressure load during operation and therefore is subjected to a possible stress 
rupture failure mode. Stress rupture is manifested as sudden, catastrophic failure of the overwrap, and is a 
function of fiber stress ratio in the composite and time at pressure. 

Stress Rupture Phenomenon 

As mentioned before, failure due to stress rupture is one of the major concerns with the usage of 
COPVs. A thorough discussion of various approaches to compute reliability under stress rupture failure is 
given in the NASA Engineering Safety Center (NESC) Carbon ITA report (Ref. 3). In recent publications, 
both Thomas and Robinson (Refs. 4 and 5) have provided models for determinations of reliability which 
included their interpretation of model parameters based on carbon/epoxy stress rupture data available at 
the time of their publications. The industry practice is to utilize these early models. The Thomas model 
was chosen to develop the stress rupture requirements in the COPV industry standard, AIAA SO-81- 
2000 (Ref. 6). 
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Thomas and Robinson models used in the COPV industry are based on Weibull distributions for 
lifetime, yet they are fundamentally different in structure. Both were generated to allow the fitting of 
stress rupture data without physics based interpretation to other composite behavior, such as strength, and 
neither is based on the fundamental fiber behavior such as single fiber failure mechanisms. The Robinson 
model is based on a Weibull power law and the Thomas model is derived from a Weibull exponential 
framework. These two frameworks were originally suggested by Reference 7. In addition, the work of 
Phoenix (Ref. 8) suggests the power law framework has more flexibility and is self-consistent with other, 
but related, composite properties such as strength. A history of these models with detailed comparison of 
the exponential and power law forms is given in the volume II of the NESC Carbon Report in Appendix L 
(Ref. 8). Due to the limitations in the Thomas and Robinson models, the Phoenix power law framework 
model was used in the NESC COPV Independent Investigation. The results from this model are 
consistent with the Robinson model as well. The remainder of the current paper focuses on the Phoenix 
classic model, and its limitations due to lack of physics and the subsequent improvements proposed by 
Phoenix which is named the fiber breakage model. A thorough discussion of these models is given in the 
Appendix L of the NESC Carbon Report Volume II. 

In the current paper the classic model of Phoenix is reviewed and its limitations are illustrated when 
proof testing of COPVs is included in the model. This is followed by a theoretical development and 
applications of a more sophisticated fiber breakage model that incorporates the proof test induced damage 
for the reliability computations of carbon COPVs. Proof testing COPVs after fabrication is conceptually 
viewed as a process of weeding out weak vessels. Vessels are most often removed because of liner 
leakage, but rarely from catastrophic failure of the composite overwrap. Nevertheless, unlike the situation 
with metal vessels, considerable damage is done by a proof test in terms of breaking carbon fibers and 
possibly strands. Classical reliability models are typically not based on fiber breakage processes in a 
way that captures the above drawbacks of proof testing. In fact, when applied to carbon/epoxy COPVs, 
the power law-Weibull version of the classical model predicts that the higher the proof test level, the 
greater the benefit in terms of improved reliability. The fiber breakage model explicitly accounts for 
the micromechanics based statistical failure processes in composites consisting of carbon fibers in 
epoxy matrix. 


Classic Model: A Review 


The usual reliability model for failure of a composite is a parametric model in a Weibull power-law 
framework that embodies a memory integral for past fiber history, < 7 ( 5 ), 0 <s<t . The distribution 
function for overall composite failure is given by 


F(t;a(°)) = l-exp^- 
\ 



t> 0 


( 1 ) 


where p is the power law breakdown exponent reflecting sensitivity of lifetime to fiber stress level, P is a 
Weibull exponent that happens to be the Weibull shape parameter for composite lifetime under fixed fiber 
stress, and a re f and t K f are characteristic stress and time scale constants for the composite material. Thus a 
stress history for a pressure vessel in stress rupture may be idealized as 


CT 


p 



0 <s<t„ 


t p <t 


( 2 ) 
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where a p is the proof level and t p is the proof hold time, typically a few minutes, and t is overall time that 
may stretch to thousands of hours. Then the distribution function for lifetime (including failure times 
during proof) is 


F (* ;a p(°)) = 1-ex P) - 
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The conditional distribution for lifetime given survival of the proof test is 
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For times much longer than the proof time (a few minutes) and well within a proof benefit timescale, t p h , 
whereby t-t p ~ t applies, this can be rewritten as 


where 
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( 6 ) 


is the timescale over which the formula applies and a proof benefit exists, since the quantity after the 
square parentheses is of order unity or smaller (typically p « 1 for carbon fiber composites). For very 
long times t p b « t the probability of failure follows 



and the proof benefit fades. 

These results can be written in terms of the conditional reliability given by o,a p ) = 1 - F[t 1 0 , 0 ^, ) . 
For times within the proof benefit range up to time t pb we have 
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Thus, according to the classic model, the reliability benefit of the proof test and how long it lasts are 
strongly related to the values of p and p , and these vary widely among fiber types. For instance, for 
Kevlar/epoxy composites, typically 1 < p < 2 and p < 30 but for carbon/epoxy composites typically 
P < 0.3 and p > 70 . 

Finally we also consider a proof test followed by a long survival time, t s , and calculate the conditional 
reliability for future lifetime of the same order of magnitude as that already survived. The probability of 
failure by time t>t s (including failures during proof and up to time t s ) 


F (* ;cT />(°)) = 1-ex Pj - 
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Noting that t p «t s <t and using the above methods, the conditional reliability given survival to time t s 
is found to be approximately 
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(ii) 


These results will be compared to those from the stochastic fiber breakage model derived below. 


Deficiencies of Classic Model 

To illustrate the effects of proof test on reliability, the relevant equations mentioned above are 
shown plotted for both Kevlar and carbon in Figure 1 . For Kevlar fibers the parameters chosen are 
p = 24, p = 1 .625 and t ref = 1 .43. These are typical values for Kevlar and are currently used in reliability 
calculations of the Orbiter COPVs. A comprehensive database of stress rupture life data for Carbon 
fibers, however, is not available. Limited data for strands and vessels are available in literature. In the 
current computations the chosen parameters are p = 72, p = 0.35 and t rq f= 0.001 . These are the parameter 
values for a typical Carbon COPV based on the limited available data. Proof test is assumed for 5 min at a 
fiber stress ratio of 0.667 and the operating stress ratio is assumed to be 0.5. The dashed lines show the 
base line case of no proof test. The dot-dashed lines show the case of effect of proof test on reliability for 
both Kevlar and carbon composites stress rupture performance. The solid lines show the conditional 
reliability where credit is taken for surviving the proof test. For Kevlar composites the proof test appears 
to have a minimal impact on the overall reliability for a period of about one year beyond which all three 
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Figure 1 . — Effect of proof test on reliability of Kevlar and carbon COPVs. 

cases are virtually indistinguishable. However, for carbon proof test appears to have substantial impact on 
the reliability for extremely long periods of time. The effects for this particular case appear to last in 
excess of 100,000 yr. Furthermore, the conditional reliability where credit is taken for surviving the proof 
test itself appears to increase by several orders of magnitude compared to the straight calculated 
reliability. Also, the proof test itself (dot-dashed line) appears to cost a significant amount of reliability in 
number of nines although it is only for a short duration. In general if one were to put a three-nine 
reliability as a minimum requirement, then carbon COPVs appear to last several times longer under stress 
than Kevlar COPVs. It is for this reason, carbon overwrapped pressure vessels are becoming more 
popular alternatives to the Kevlar counter parts in space applications. 

In conclusion, if the classic model is correct, the proof test benefit for Kevlar fiber is relatively short 
lived while that for the carbon fiber could last for any conceivable lifetime of the vessel, provided that 
p > 70 . However, the model is based on purely curve fitting to statistical stress rupture data. Any damage 
that may have resulted during proof test is not accounted for in the model. Furthermore the model in its 
current form indicates that higher the proof stress ratio, the better it is for reliability provided the vessels 
survive proof test. As the proof test stress ratio increases composites tend to accumulate damage in terms 
of fiber breakage which will have adverse effects to its performance. By considering the process of fiber 
failure more explicitly, we investigate whether the classic model is a realistic reflection of the failure 
process in the case of carbon/epoxy composites or whether it is a misleading idealization. 


Fiber Breakage Model Development 

The fiber breakage model theoretical development is discussed in detail in References 8 and 9. For 
sake of completeness, this discussion is included in Appendix A. Key equations that are needed for 
discussion here are shown below while the details are deferred to the Appendix A. The life time 
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distribution function (probability of failure for T< t for composite under stress a which has gone through 
a proof test at g p a hold period of t p is given by Equation (A.3 1 ) in Appendix A: 



where 6 V composite reference strength (same as a ref ), k is critical fiber break cluster size for 
catastrophic failure, k p is fiber break cluster size after proof, q is Weibull shape parameter for fiber 
strength, and 0 is the matrix creep exponent, a is given by the product kq . The fiber break cluster size 
after proof k p is given by Equations (A.29.a) and (A.29.b) in Appendix A. 



(13) 


The conditional reliability given that any weak vessels removed by proof test is given by the following 
equation (App. A, Eq. (A.35)) 



Finally for the case where the vessel has not undergone any proof test the reliability is given by 


. i \<* ( /i . 777 . \ e ^ 


rt(f|CT)»expj-(a/CT K ) a Jl + (f/f c ) >, 0<t 


Equations (5) to (7) can now be used to perform comparison studies to classic model for COPVs. 


Comparison of Classic Model Versus the Fiber Breakage Model 

To compare the predictions of classic model to fiber breakage model, a T1000G COPV with a 
nominal stress ratio of 50 percent is considered. The proof factor (defined as the ratio of fiber stress at 
proof to fiber stress at operating conditions) is 1 .5 which gives a stress ratio of 0.667. These are typical 
values for such vessels. The model parameters for T1000G fiber are taken as: 


NASA/TM— 201 0-21 583 1 


7 



Classic model: 


p = 114;* cref =0.001hr;P = 0.22; (16) 

Fiber breakage model: 

c = 5 \t c = 0.01hr;£ = 5 \t p = 0.0833 hr;d = k£ (17) 

Predictions of reliability of T1000G COPV expressed as “number of nines” are shown plotted as a 
function of time for both the models in Figure 2. Calculations are shown for both no proof test case and 
the case where survival of a proof test conducted at 66.7 percent stress ratio for 5 min. 

Both models predict almost identical reliability when proof is not considered, but the classic model 
predicts vastly improved reliability when conditional survival of proof test is taken into account. 

However, the fiber breakage model predicts much lower reliabilities when proof test effects are included 
and the reliabilities are based upon conditional survival of proof test. A loss of two nines in reliability is 
possible due to proof test as indicated in Figure 2. The reason for this divergent behavior is due to the fact 
that the fiber breakage model takes into account fiber damage that takes place during the proof test. It is 
for this reason recommended that the level of proof test be carefully determined such that excessive 
damage to fibers is avoided. 

Parameter Sensitivities and Uncertainties 

The sensitivity of reliability to changes in various parameters is shown in Figure 3. Here each 
parameter is varied with respect to its nominal value. The nominal values chosen for the parameters are 
for the T1000G carbon fiber and are given by Eq. (10). 

£ = 5 \t c = 0.01hr;£ = 5 ;t p = 0.0833hr;A = *£ ; 



Time in Hours 

Figure 2. — Comparison of classic and fiber breakage model reliability predictions 
for a T1000G COPV operating at a stress ratio of 50 percent. 
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Stress ratio considered is 50 percent and proof stress ratio is 75 percent. Each parameter is made 
dimensionless by dividing with nominal values and the dimensionless parameters are varied from 0.7 to 
1 .3, one at a time. The reliability of the vessels conditional on proof test survival is computed for 10,000 
hr and is shown plotted versus all the parameters. 

It is clear from the figure that the Reliability shows a high degree of sensitivity to changes in the 
parameters £ followed by k ’ and 0. The parameters t c and t p have very little or no effect on the reliability. 
Determination of the three most important parameters to a high degree of accuracy is therefore essential 
to assess the reliability of the vessel accurately. 

In general when dealing with reliability one must assess and account for two types of uncertainties: 
epistemic uncertainty and aleatory uncertainty. Epistemic uncertainty arises due to lack of knowledge, or 
insufficient data. Such uncertainties can be reduced over time with more experiments and experience. On 
the other hand aleatory uncertainty is the so called physical variability that is present in almost all the 
systems and is unlikely to be reducible. Various publications address these two types of uncertainty and 
the importance of separating and resolving them (Refs. 9 and 10). 

These uncertainties are especially critical for the COPVs because of the risks and costs involved. 
Parameter or epistemic uncertainty can be reduced by collecting more data and by developing better 
analytical models. Here we attempt to show how the model form uncertainties affect the reliability 
estimates via Monte Carlo simulations of the conditional probability of survival. As an example a typical 
vessel made of T1000G is chosen which has been successfully on flight for the past 7 yr. The parameters 
that are considered uncertain with the respective distribution statistics are given Table I. These values are 
taken from the References 3 and 1 1 . An observation of interest is the fact that the variability in parameters 
is rather high. This is due to very limited stress rupture data that is available for this class of materials 
(Ref. 1 1). The variability in parameters is established based on this limited data and one can expect to 
reduce this when more data is available thus improving the reliability considerably. 

The proof stress ratio is taken to be 0.75 and proof test time t p = 5 min. The stress ratio for the vessel 
is taken as 0.478. These values are typical of the highest loaded COPV on the ISS. Reliability conditional 
upon surviving successfully the past history of 7 yr is computed using both classic model and the fiber 
breakage model and are reported in Table II. The calculations are repeated for missions up to the year 
2016, 2020 and 2030 starting from year 2009 which represents a past history of 7 yr. In these calculations 
it has been assumed that the vessel is at the maximum expected operating pressure (MEOP). 



Parameters (£,k , ,t c ,t p ,6) 

Figure 3. — Reliability of a T1000G Carbon fiber COPVs for 10,000 hr as 
a function of various model parameters. 
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However, typical COP Vs experience a ramp reduction in pressure with time as the gas is used. The 
assumption that the COPVs maintain MEOP is conservative for operation beyond the conditional point 
and unconservative for the past history. However, in this example case, the assumption of MEOP 
operation before the conditional point did not contribute significantly to the reliability and was neglected. 


TABLE I.— PARAMETERS USED FOR CLASSIC AND FIBER BREAKAGE 
MODELS IN THE UNCERTAINTY SIMULATION STUDIES 


| Classic Model Parameters | 

| Parameter 

Mean 

cov 

Distribution 


114 

0.25 

Lognormal 

p 

0.22 

0.3 

Lognormal 

^ref 

0.001 

1.0 

Lognormal 

| Fiber Breakage Model Parameters | 

Parameter 

Mean 

COV 

Distribution 

9 

5 

0.25 

Lognormal 

0 

0.11 

0.25 

Lognormal 

£ 

5 

0.3 

Lognormal 

t c 

0.01 

1.0 

Lognormal 


TABLE IL— RELIABILITY PREDICTIONS FOR A TYPICAL CARBON COPV ON SERVICE 
CONDITIONAL ON SURVIVAL OF 7 YR OF PAST HISTORY AT MEOP TO OUT YEARS 


2500 psi operating 
(Stress ratio = 0.478) 

2016 

(70,080 hr) 

2020 

(105,120 hr) 

2030 

(192,720 hr) 

Point 

0.999943 

0.999925 

0.99988 

Mean 

0.999521 

0.999279 

0.99894 

95% 

0.99768 

0.99738 

0.99615 


Three different estimates point, mean and the 95 percent confidence are computed by running 
thousands of Monte Carlo simulations Reliabilities are conditional on successful past history of 7 yr. As 
seen in the table, the classic model does predict a slightly higher reliability (three effects. The classic 
model predictions are based on neglecting the proof test effects because if proof test were included, the 
reliability would be very high as suggested in Figure 2. The 95 percent confidence estimates, however, 
fall slightly short off the three nine mark in this example. As can be expected the reliabilities can be 
increased by reducing the pressure if the three nine reliability to the end of the mission is absolutely 
critical to the program. 

The second example in this section pertains to hypothetical reliability assessment of next generation 
launch vehicle. Here the various models and predictions are assessed and compared for a system of 
COPVs made of T1000G fiber. Because the vessels have not been manufactured, the exact stress ratio is 
not known. For purposes of this initial analysis, a proof factor of 1 .25 and a design burst factor (DBF) of 
2.0 are assumed. The stress ratio is given approximately by. In Figure 4, the DBF is varied from 1.5 to 3.0 
and proof factor is scaled appropriately from the 2.0 DBF case. Reliability is computed for both classic 
and fiber breakage models. Tables III and IV show the details of the classic and fiber breakage model 
parameter values chosen in the reliability computations. 


TABLE m.— CLASSIC MODEL CASE 1 AND 2 PARAMETER VALUES 


! Classic Model 1 | 

Parameter 


Nominal value 

P 


114 

p 


0.22 

^ref* 


0.001 

Classic Model 2 

Parameter 


Nominal value 

P 


72 

p 


0.35 

^ref- 


0.001 


NASA/TM— 201 0-21 583 1 


10 
















TABLE IV.— FIBER BREAKAGE MODEL PARAMETER VALUES 


Fiber Breakage Model Parameters 
(High variability) 

S 

5 

0.25 

Lognormal 

0 

0.11 

0.25 

Lognormal 

k 

5 

0.25 

Lognormal 

tc 

0.01 

0.3 

Lognormal 


Fiber Breakage Model Parameters 



(Low variability) 


S 

5 

0.05 

Lognormal 

0 

0.11 

0.05 

Lognormal 

k 

5 

0.05 

Lognormal 

t c 

0.01 

0.05 

Lognormal 


The simulation results based on the above four models are shown in Figure 4. The two classic models 
are chosen to represent the case of no prior proof test and the case where proof stress induced damage is 
accounted. Thus the classic model case 2 parameters are chosen to simulate reduced performance due to 
possible proof test fiber damage. Two fiber breakage models one with a high degree of variability and one 
with low variability are considered to represent the current state of knowledge and a future state of 
knowledge with better understanding of carbon stress rupture behavior and the availability of a 
comprehensive stress rupture database. Thus, the high variability of fiber breakage model represents the 
current status of the limited data and knowledge of carbon stress rupture behavior. The low variability 
case represents what we can expect for reliability in future when our knowledge base of stress rupture 
behavior for T1000G fiber is enhanced through substantial stress rupture testing and establishing a 
comprehensive data base. As can be seen from the figure, for a COPV system of 13 vessels, and a DBF of 
2, the reliability can vary from a little over one nine to five nines depending upon the model chosen. The 
uncertainty in reliability high lights the overwhelming need for a formal development of an experimental 
stress rupture testing program 



Figure 4. — Various model predictions of reliability of a system of vessels on board a launch 
vehicles for a total mission time of 1 yr. 
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Summary 

Stress rupture life and reliability of COPVs for space applications are assessed by using two different 
models with some illustrative examples. The classic model while simple to use has the short comings of 
not able to represent directly the proof test induced fiber damage. The conditional reliability based upon 
survival of proof test shows unrealistically high value when classic model is used. The reliability appears 
to increase as the proof stress ratio increases and its effect lasts for thousands of years. This is clearly an 
artifact of the model which is purely phenomenological model and lacks physics. Proof tests performed at 
higher loads induce fiber damage due to scatter in fiber strength and that is not captured in this 
phenomenological model. In order to address this micromechanics based fiber breakage model is 
proposed. The model accounts for possible fiber damage due to proof test. However, this model gives 
overly conservative reliability estimates. Some of this conservatism is due to lack of sufficient data for the 
T1000G carbon fiber that is considered in the simulations. This introduces huge uncertainty in the 
parameter estimates which translate into substantial reductions in reliability estimates. 

The wide range of reliability estimates depending upon which model is used warrants immediate 
attention to the development of a formal experimental test plan to reduce this uncertainty and narrow 
down the differences between the models or formulate better analytical models based on the 
comprehensive test data. Furthermore, the fiber breakage model is based on progressive damage theory 
which has not been rigorously tested and validated. Currently the NESC is planning to address some of 
these issues by a detailed stress rupture testing that include proofed strands and vessels as well as un- 
proofed test articles. It is anticipated that the tests will provide more clarity to various models and help 
reduce the uncertainties as well as improve some of the models. The results will be reported in subsequent 
papers. 
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Appendix A 

Stochastic Fiber Breakage Model Development 

Micromechanics and Statistics of Fiber Failure Process 

To gain a deeper understanding of the issue we first discuss the failure process in the composite in 
terms of fiber stress redistribution around fiber breaks (i.e., load sharing) and local break cluster 
formation among a large number of fiber elements of length, 8 e , which is the characteristic elastic load 
transfer length from a broken fiber to an intact neighbor. This length will be compared later to a length 
5(t) that grows with time beyond a characteristic time scale t c for matrix creep in shear around initial and 
evolving fiber breaks in a cluster. Note that t c is of the same order of magnitude as the “proof hold time”, 
t„, in a proof test. Also 8 e depends on many mechanical and geometric quantities, the most important 
being the fiber diameter, df, the fiber Young’s modulus, Ef, the matrix shear modulus, G m , and the 
effective width of matrix b etween the fiber s. w m , which depends on the fiber volume fraction, Vf. For fully 
elastic behavior 8 e ~ 2d f jEyw^ j[G m d^ and since 140< / G m < 200, w m jd f » 0.3 , and 

d f « 5 pm we obtain 8 e » 20 d f » 0.1 mm . 

There are a large number of such fiber elements in the overwrap of a typical carbon/epoxy pressure 
vessel. A single yam or tow typically has 12,000 fibers, and one wrap covers about 2.5 mm of width 
around the vessel circumference and is 0.25 mm thick. For a COPV of diameter 75 cm, the overwrap 
thickness may be 4 mm corresponding to 8 overwrap layers. Thus for a typical carbon/epoxy pressure 
vessel there are of the order of 10 12 fiber elements of length S e in the overwrap. Later in the analysis, the 
volume, V, of the overwrap will be the number of such fiber elements of length S e that it contains. 

The probability that a fiber element of length 8 e breaks under stress a is well approximated by a 
Weibull distribution of the form 


F s e ( CT ) = 1 - expj-(8 e // 0 )(ct/ct /o f j, 

= 1 — expj — (cr/cr s< n, ct>0 


(A.l) 


where l 0 is a typical tension test gage length for fiber testing (e.g., 1 cm) with associated Weibull scale 
parameter cr /o , and shape parameter £ estimated from such testing. Also a s ^ is the Weibull scale 
parameter for the strength of a fiber element of length, 8 e . Implicit in the above formula is the scaling of 
fiber strength with length following 


=<*/ 0 (8 e A)) ( A - 2 ) 

which reflects the fact that that carbon fibers typically follow Weibull weakest-link statistics. Typically 
C, » 5 for carbon fibers so the length effect is fairly strong. A ten-fold decrease in length increases the 
strength by about 58 percent. A useful approximation for (a), especially in the lower tail, is 

F s e (°) ~ (°/ CT s e f > 0 < a < a 5e (A.3) 

The characteristic strength, cr 8 , is typically about 3 times the effective stress level, a v , or a s « 3oy . 

Next we consider the fiber stress level a « ct 8 and proceed towards calculating the probability of 
formation of a cluster of size k emanating from a given fiber element of length 8 e in the composite. Note 
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that the number of such “initial” breaks is a random variable that is binomially distributed with 
parameters V and (ct). The mean number of such breaks is h(ct) « VF d ^ (ct), and thus, does not 
depend on the calculated value of 5 e since a « ct 8 so that breaks are far apart. In fact only a small 
proportion of the fiber elements fail under the applied stress, ct. For instance, for ct/ct f = 0.50 , or a fiber 
stress ratio of 50 percent, we have ct 8 » 3a v ~ 6a and for C, » 5 the probability of failure of a given fiber 
element is approximately (ct/ct 5 ^ ) « ( 1/6 ) 5 » 1 .3x 10 -4 . Thus, about one in every 8,000 fiber elements 
fails under stress ct and along a filler the breaks are about 78 cm apart. While this may seem to be a large 
spacing, the number of such broken fiber elements in the overwrap is roughly h(ct) « VF S (ct) 

~ (l.3x 10 _4 ^10 12 ) = 1.3x 10 8 , still a very large number. 

All of this assumes a homogeneous or uniform stress field. If there are areas of stress concentration, 
say near a boss where failure is more likely to occur, then a relatively smaller volume of material might 
be at a higher stress concentration, and thus, the breaks may be more concentrated. There still will be a 
large number of fiber elements in this region of higher stress and the breaks will be more frequent, i.e., 
more closely spaced. By themselves these initial breaks obviously do not cause composite failure and 
some micromechanical interaction is necessary to generate break clusters that eventually become 
unstable. This is pursued next. 

Derivation of Composite Strength Distribution 

The failure process can be characterized in simplified terms as follows: First is the growth of clusters 
of fiber breaks. In order to fail the composite one of these clusters must grow to a critical size. Locally the 
formation of a cluster begins with failure of a single fiber element of length 8 e under stress ct followed by 
failure of a neighboring element under the stress concentration K x a , and another neighbor under higher 
stress concentration K 2 a (since there are now two adjacent breaks) and so on until the size reaches k. 

(The stress concentration factors grow unbounded and satisfy 1 < K x < K 2 < ... .) The probability of 
formation of such a cluster is approximately 

W k (a) « Ck F &e (a)F &e (K,a)F &e ( K 2 a )• ■ ■ F, ^ (K k _,a) 

* c k {K x K 2 • • • K k _ x f ( ct / ct 8j 

where c k is a combinatorial factor capturing all the possible configurations (in terms of sequence of fiber 
breaks) that a cluster can have. For instance, in a planar array of fibers with load shedding primarily to 
nearest surviving fibers (one on each side) it can be shown that 

KjvJl + nj/4, j = 0,1,2,- •• (A.5) 


and 


c k *2 k ~ l (A. 6) 

The factor c k = 2 k ~ x captures the fact that, except for the failure of the given fiber element to start the 
cluster, there are typically two overloaded neighbors next to the growing cluster at any growth step. At 
applied fiber stress level, ct, the critical value of k that triggers a catastrophic cascade is called and 
satisfies 

K k'-\°<°8 e < K k° ( A - 7 ) 
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since beyond this value of k, all neighboring fiber elements to the cluster will be overloaded well beyond 
their characteristic strengths . Once this occurs, their individual probabilities of failure now exceed 
0.632 and with further cluster growth rapidly approach unity so that the cascade cannot be stopped. 

In a hexagonal or ‘random’ array of fibers the ideas are similar but c* will grow much more rapidly 
than in above planar case and will involve products of increasing numbers of fibers around the periphery 
of a growing approximately circular cluster. However the values of K x , K 2 ,- • •, K k ,- • • will grow more 
slowly. For instance, it can be shown that 

Kj * yJl + D/n, j = 0,1,2,-- (A.8) 


and 


c k * ** TI(^ + 0 ( A - 9 ) 

7=1 

where D is approximately the cluster diameter, and j = nD 1 / 4 is approximately the number of fiber 
breaks in the cluster. 

It can be shown that the distribution function for the strength of the composite at a fiber stress level a 
is approximately 


which can be written as the Weibull distribution 

//p(a) « 1-exp 



(A. 10) 


(A. 11) 


where 

= °8 [Vc k )"^ (K x K 2 ■ ■ ■ K k ] y' /k (A. 12) 

and 


a. = kC, 


(A. 13) 


and where 


^_,a K <a 8e 


<K,-a 

k 


V 


(A. 14) 


which is equivalent to k = k~ . 

In the above description we ignore the fact that k, increases slowly as c decreases (whereby the 
probability of composite failure typically decreases by more than an order of magnitude with each integer 
valued increase in k,). Instead we have chosen the particular value, k , called the critical cluster size for 
catastrophic fracture, when the stress level is equal to the Weibull scale parameter for composite strength. 
In this case k satisfies 




(A. 15) 
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This equation implies that the predicted strength distribution for the composite is not actually a straight 
line but is slightly curved on a Weibull plot and becomes steeper at lower stresses implying slightly lower 
probabilities of failure than a straight Weibull line would predict. Typically the change in slope is so 
small that it cannot be detected in pressure vessel testing and carbon/epoxy strand testing with typical 
number of test samples. However it has been very clearly detected in the experimental testing of seven 
fiber micro-composites on a repeatable basis (Refs. A.l and A.2) and is seen also in Monte Carlo 
simulations (Ref. A.3). Furthermore there is some evidence that as structures become larger their 
variability in strength decreases, which is consistent with the increasing steepness feature above. 


Composite Lifetime Distribution in Stress Rupture in Absence of Initial Proof Test 

Stress rupture in carbon fiber composites arises mainly because the matrix creeps in shear around 
fiber breaks. A useful model is that this creep follows a power law creep compliance of the form 




t> 0 


(A. 16) 


where J mfi is the instantaneous elastic creep compliance (J me = \IG me ), where G m e is the instantaneous 
elastic shear modulus), t c is the characteristic time for creep to occur (at which time the compliance J m (t) 
has effectively doubled), and 0 is the creep exponent. The creep exponent is a crucial parameter and 
depends on such factors as the matrix and interface chemistry in terms of adhesion to the fiber, fiber 
volume fraction, and temperature — to name the most important influences. Typically 0.1 < 0 < 0.5 for 
epoxies, and note that as a reference point the value 0 = 1 corresponds to a Newtonian viscous material. 
We note that nonlinear forms of the creep function described in terms of visco-plasticity introduce a 
dependence on shear stress level between the fibers raised to some power larger than unity (i.e., 3 or 4). 
However this modifies the analysis in only modest ways and does not influence the conclusions. 

It can be shown that, under the power-law creep compliance, the characteristic load transfer length 
around a fiber break is time dependent following 

8(<) = fW 1+ ('/<X 

where 5 e is again the instantaneous elastic load transfer length. Thus as t increases, this overload length 
grows so that increasing numbers of flaws in fibers adjacent to a break cluster of size k are exposed to 
stress K k a and can fail despite surviving the original stress level a. 

Susceptibility to such delayed failure is greatly enhanced for carbon fibers with high variability in 
strength, i.e., low values for the Weibull shape parameter, £, for strength. Individual carbon fibers are 
virtually immune to stress rupture compared to that seen in carbon/epoxy composites. Essentially single 
fibers under constant stress fail on loading or do not fail at all. Of course, in the composite matrix creep 
makes it possible for the stresses on fibers next to broken fibers to increase and so the main cause of their 
failure is increasing stress around growing clusters. 

In stress rupture at a fixed stress level, o<6 v , and assuming no previous stress excursion above this 
level (such as a proof test), the lifetime distribution function can be derived as a modification of the case 
above. The distribution function for composite lifetime follows 

H- K (r;a)*l-exp[-F^(r;a)], t»t c (A.18) 


t » 0 


t»L 


(A.17) 
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where W~ (f;a) is a characteristic distribution function given in simplified form by 


"CtfcK. 2 -K t _ x 


i) ; ( a / a s. r(V i+ ( # Ac) 0 ) 

n*Kmr». 


(A.19) 


This form reflects the fact that at a given element location the creep-rupture process requires 
failure of that element simply due to the applied stress, followed by the failure of some element out of 
2 Jl + [t/t c ) neighboring fiber elements along the adjacent fibers exposed to the overload Kf > , where 
the number exposed grows over time as the load transfer lengths grows. This process continues as more 
and more elements fail in time and ultimately these time dependent failures in a cluster number k - 1 . 
Fiber breaks occur sequentially so there is some time lag for growth of the new overload length at each 
new fiber failure site, but this is not reflected in the above formula where time is just the original time, t. 
Simulations show, however, that the effect is small compared to the long timescales involved (Ref. A.3). 
A more refined analysis would lead to the need for the multiplicative factor of the form, 


®(*)=r(9/2+i)/r (it- 1 ) 0 / 2+ 1 


(A.20) 


where T(z) is the usual gamma function. However, typically 0.05 < 0 < 0.25 , and unless k is extremely 
large, this factor is not significantly different from unity and can be neglected. 

The resulting Weibull approximation is 


H r [t;c)» 1-exp 



t»t c 


(A.21) 


where 


p= (it -l)0/2 


and 


W - *c( ct /^k) P 


are the respective Weibull shape and scale parameters for lifetime. Also 


P = 



(A.22) 


(A.23) 


(A.24) 


and 


p = 2 C/0 (A.25) 

are the power-law exponents for lifetime versus stress level, being the effective value and large cluster 
limit, respectively. 

In summary, in a stress rupture setting the applied stress, a, is much smaller than the scale parameter 
for tensile strength, a v , so that the probability of failure on initial loading, H v ( cr) , is easily satisfies the 
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reliability requirement. However, probabilities of failure at much larger times l » l c become the key 
concern. Essentially, the first failure at some location (necessary to start the matrix creep process) is due 
to a flaw in the fiber element of elastic length 5 = 5 e , and the remaining failures occur due to the time 
dependent matrix creep process since the overload lengths quickly grow beyond 5 e for t»t c . Since we 
assume no initial overload (such as a proof test), the initial clusters of breaks due to the initial loading are 
automatically consumed by those that occur in time since any fiber that fails initially is automatically 
included in time dependent failures. Another important fact inherited from linear viscoelasticity is that k 
for stress rupture is almost the same as k for strength at times near zero, so we do not distinguish 
between the two. 


Effect of Proof Testing 

As mentioned, proof tests are often applied to pressure vessels to some fiber stress level, o p , in order 
to filter out weak vessels. For all metal pressure vessels this process can be argued to be all beneficial 
with no drawbacks, i.e., no new damage is introduced; however the situation is far less clear in composite 
over-wrapped pressure vessels. The difference is that, because of the proof test, many fiber elements will 
now fail under the higher stress o^, that would not have failed under the lower stress-rupture stress level, 
a. These additional failures now provide many additional locations for subsequent cluster growth in 
stress-rupture. 

For instance, the number n(a p ) of fiber breaks at the proof stress level o p divided by the number n(a) 
at the stress-rupture level a is given by 

n { c p)/ n (°) = ( G p/ 6 vf =( CJ P / CT f ( A - 26 ) 

For the typical value C, = 5 and = 1 .5 we obtain «(1.5 ct)/«(ct) = (l.5) 5 = 7.6 . Thus, there are 
7.6 times as many single fiber breaks or ‘singlets’ due to the proof test as without the proof test and these 
additional singlets now become many more seeds for stress-rupture compared to those that would occur 
without the proof test. The situation is made worse, however, since some singlets among the expanded 
group (from the proof test) will generate many failing neighbors due to the overload K x o p that would not 
occur due to the milder overload Kfs , thus forming additional ‘doublets’ of failed fibers. Furthermore 
some of the neighbors to these new doublets will fail under K 2 o p but not under K 2 o , thus forming 
additional ‘triplets’, and so on. In essence, as a result of the proof test, there will be a distribution of fiber 
break multiplets of various sizes, and the numbers of such multiplets of each size will be much larger than 
would have occurred without the proof test. 

A proof test to stress level, o p , will eliminate any vessel that develops a cluster that reaches critical 
size, which for simplicity we take to be k , i.e., a weak vessel. Unfortunately it suddenly creates many 
clusters that would not occur otherwise. This is very different from the model for metals or ceramics, 
where the proof test may eliminate any vessels with flaws or cracks larger than critical, but anything 
smaller will actually be beneficially blunted and no new ones will be created. In ceramics the model is 
similar, though each smaller crack may grow slightly. It must be noted, however, that the fiber itself 
benefits from the proof test in that the vast majority of the fiber elements that survive the proof stress 
level, o p , cannot fail until the stress level in the future exceeds o p due to some overload Kp next to a 
sufficiently large cluster. 

Next we consider an initial proof test to fiber stress level, o p , over time 0 < t < t (typically involving 
a ‘proof hold time’ of a few minutes) and let W~ (t;o,o p j be the characteristic distribution function for 
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lifetime extending the concepts developed earlier, except we now have two stress levels, a and a p . We 
then have 


H r (t;a) * 1 -exp{ -VW- ( t;a p ,a p )} , 0 <t<t p 


(A.27) 


the time period up to t p , where 


w t, P (c®.®,)” ( k < k > ■ ■ ■ k >-< t kk r { + ('/<c) 9 1 ' - o<,<, r < a - 28 > 


This is identical to the form earlier except that now <j = a p . 

We now consider the effect of the stress decreasing from a p to o for t > t p . Although further growth 
of the overload length on fibers next to a break cluster continues at approximately the same rate, the 
overstress itself is now reduced by the factor o p . Thus the rate at which new breakable flaws are 
encountered is reduced nominally through the Weibull factor (a/a p ) . The overall effect is nominally 
equivalent to the stress remaining at level a p , but the rate of increase of the overload length being slowed 
in time by the factor (a/a p ) . This characterization, however, is not entirely correct since all fiber 
elements have been proof tested to a p and any further fiber element failure requires a growing overload 
length around some cluster of size j, such that Kp > a p . The critical size for such a cluster is denoted 
k p = k p (o/<j p ) and satisfies 


K k p -P<v p < K kV or K k p -i<v p /v< K k p 


(A.29.a) 


Thus any cluster of size k p after the proof test will continue to grow, and smaller ones will not. We see 
that cr p /cj > 1 requires k p > 1 and using K k « sj\ + nk/4, k = 1,2,3... we can see that 


( CT „ Y ] 


v a y 


-l 


+ i 

y y 


(A.29.b) 


where [z] refers to the integer part of z. In computations, it suffices simply to use 


k P K 


n 


( 'f ] 


v y 


-1 


+ 1 


(A.29.C) 


Thus for t p < t we obtain 


W kr (,^P H (k,K 2 . . . K u f (a,/a 5> f [ fpp' 




^ 1 + (^Ac) + ( <J / <J pf('] 1 + ( t / t cf -^j 1 + ( t p/ t c) 


1 

Jj 


k-k„ 


(A.30) 


Using again (A.23) for a v , we reduce (A.41) to 
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(A.31) 


K /*V T \J 1 + Vp Ac f y 
*\L + { a / a pfij 1 +(*/*€?, 


W k ~ k P 


, t>t D 


and finally 


H v ) * 1 - ex pj|( CT P / d F )“ ^ 1 + (^Ac)' 


\ k ~' 


1 + ( CT / CT pf^( 1 + (^c) e )/(l + (^Ac) e )-l| || . l>t p 


!*-*/> ) 


(A.32) 


Conditional Reliability After a Proof Test 

Next we must evaluate the conditional reliability given that any weak vessels removed by the proof 
test. The conditional reliability given survival, is 


R(t I a * a p) = ex P {-^ \w itp (t,o,o p ) - W k p (t p ;o p ,a p )] 
where the characteristic distribution function in time, t, beyond the proof time is 


(A.33) 


^ip (^p ) * ( a p/ 6 r T {I 1 + (*pA f y 

x | 1+ ( a / CT pf -■ 


w k ~ k p 


(A.34.a) 


. t p <* 


whereas 


) * ( a P / d v f [ ^ 1 + (^Ac)' 


\* -1 


0 <t<t„ 


(A.34.b) 


Thus the conditional reliability, (A.44), becomes 

*(* 1 ct,ct p ) » exp( -(ct p /6 V )“ f. \jl+(t p /t c ) 




l+ ( a / a pf ^( 1+ ( f Ac) e )/( 1+ (^Ac) ) _1 


l 1 

A 


k-k„ \] 


(A.35) 


-1 


y 


K o <t p <t 


Two particular cases are of interest: If the proof time t p is such that ( t p j t c j « 1 , then we approximately 
have 
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[ 




i?(f|CT,CT p )»expj- 

ll 


K a v j 


On the other hand if ( t p jt c ) « 1 , then approximately 


t kp 


-1 


ll 

K 0 <t <t (A.36) 
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*(*l CT » CT i»)" e *P?- 
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v ct k y 


>(M / 2 


1 + 


y >£ 


v CT P y 


V2' 


1 + 


V 9 

\ l c J 


\T kp 


-i 


-i 


il 

K t c <t 

U 


(A.37) 


We can see that the behavior of the reliability in (A.48) depends on the magnitude of the time t relative to 
a transition time t p b that solves 


( \ 

p 

/ 



0 | 

a 


1 

— =A 

|l + 


-1 

< c p J 


V2]j 

V 


U J 

J 


= 1 


(A.38) 


or 


1/0 


tp* =t \ 2(a p /cf +2j2(<y p /cyf 
Then for moderate times such that t c « t « t p b we can expand the (A.48) to obtain 

nT /2 ) 


(A.39) 


R (t I v,v p ) * expj -2 i/2 1 (£ - k p ) ^ 


T 


T 

v J 


r v 


v CT *y 


r . \i 


\}c 


}, t c «t«i p>b 


(A.40) 


'J 'J 


where again p = 2^/0 . On the other hand for extremely long times beyond the threshold time t p b , we 
obtain from (A.48) a different expansion of the form 


^(t|a,CT /) ) = exp-{-2 

l 


>*/2-l 


f ~ 


V a r J 


f V 


v G />y 




t c Jj 


,jl 

>, i P jb«* 

J 


(A.41) 


where 


p p =(k-k p )e /2 

Proof Test and Survival to Some Long Time 


(A.42) 


Finally we consider the conditional reliability for the case where the there was originally a proof test 
but the pressure vessel has also survived for a long time t s »t c . Then the reliability can be shown to be 
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(A.43) 


For (t p jt c ^ » 1 we can simplify the above to 


^(tla.CTp^^JwexpJ- 
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( ^ A a 


V ct k y 
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(A.44) 


Also for intermediate times satisfying t c « t « t b this expands to 


R(t | c,c p ) » exp^ -2 k/l ~' (£ - k p ^ 


\“ 


p_ 

v ) 


f Y 


y c p j 


\ l c 


<1 
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i? , a 0 / 2 m 
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J-, t c «t«t 
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(A.45) 


However, for extremely long times i p h « t we have 


I CT >°/> ) = expj -2 k ! 2 1 


\“ 


k CT r y 


y y 


v c ^ y 


aT 


/c Jj 


vA y 


\1 


-1 


K ^ «t 
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(A.46) 


where 


$ p ={k-k p yi2 


(A.47) 


Comparison of Stochastic Fiber Breakage Model and 
Classic Model in Case of a Proof Test 

We can compare the results of the classic model and the stochastic fiber breakage model to see what 
the key differences are. For the classic model we assume t p » t ref « t c and naturally let cr ref = a v and 
a = pp. Then 
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Note that the expression (A.59) for the conditional reliability contains the fundamental factor 


<I> p =P 
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In the classical model <£> p is linear in time, t, so early on when, 0 < t p «t« t pb , we find that O p is 
very small and captures the reliability improvement due to the proof test. 

In the stochastic fiber breakage model just developed and with a short proof test hold, 0 <t « t c 
the situation is different in that we have 


R(t | c,c p ) » expjj -(*- k p 
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Thus the factor ® p now becomes a very different factor, d> p , given by 
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(A.52) 


Unfortunately 0 is a very small number (e.g., 0.08 < 0 < 0.16 ) for carbon fiber/epoxy composites. Thus 
the benefit of the proof test is very short-lived as d> p rapidly approaches unity. Then the reliability 
rapidly becomes that of surviving the proof test itself. 

We also notice that the difference between the stochastic fiber breakage model and the classic model 
tends to disappear for larger values of 0. In fact, for 0 = 2 and noting p » a/p « ^/(2^/0) « kd/l = k 
we see that <t> p « <t> p so the two results become the same. Typically 0 growing larger corresponds to p 
becoming smaller, i.e. the composite is more susceptible to stress rupture overall. However the 
differences are most marked for the case of large p. 
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